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PREFACE 


The  essential  ideas  of  this  paper,  the  introduction  of  mixed 
strategies  and  the  modification  of  Everett's  [5]  and  Pugh's  [9] 
double  Lagrange  multiplier  method  for  constrained  min-max  problems, 
have  had  a  checkered  career.,  .  . 

These  ideas  arose  from  a  computational  device  that  completed 
the  solution  (or  appeared  to  complete  the  solution)  of  the  ABM/ 
Shelter  deployment  problem  formulated  by  Robert  Kupperman.  The 
device  was  a  powerful  one  and  I  promptly  applied  it  to  a  much  more 
complex  and  interesting  version  of  the  same  problem.  The  resulting 
offense  and  defense  solutions  were  thoroughly  convincing. 

What  has  followed  this  initial  discovery  and  application  has 
been  a  succession  of  attempts  to  explain  why  the  device,  now 
christened  the  modified  double  Lagrange  multiplier  method,  apparently 
works.  The  first  attempt  appears  to  me  now  to  have  been  quite 
naive.  It  foundered  with  the  discovery  that  a  presumed  minimum  was, 
in  fact,  a  local  maximum.  The  second  trip  to  the  drawing  board 
produced  a  saddle-point  theorem  which  works  if  the  number  of  feasible 
strategies  for  the  maximizing  player  is  finite.  However,  all  of 
my  attempts  to  extend  this  theorem  to  a  continuum  of  strategies 
exploded  in  one  way  or  another.  Worse,  the  saddle-point  theorem 
could  not  provide  a  basis  for  my  method  unless  it  could  be  so 
extended. 

This  paper  follows  an  entirely  different  route  from  these 
earlier  attempts,  I  am  frankly  indebted  to  Owen  [  J  whose  paper 
contains  the  essential  step  I  have  borroweu  and  put  to  my  own  use. 

The  results  that  I  have  (I  hope)  proved  arc  not  as  strong  as  I  had 
initially  expected.  Nevertheless  they  should  provide  a  sufficient 


ii 


basis  for  the  more  grandiose  applications  of  the  modified  double 
Lagrange  multiplier  method  that  I  have  attempted  elsewhere. 

I  am  greatly  indebted  to  a  number  of  my  colleagues  who 
patiently  studied  the  confused  logic  and  opaque  language  of  several 
previous  drafts  of  this  paper.  In  particular,  the  comments  of 
James  T.  McGill  of  IDA  and  James  E.  Falk  of  RAC  led  to  the  correction 
of  several  errors  and  helped  me  to  improve  the  exposition  throughout . 
The  research  reported  in  this  paper  was  supported  by  the  Office  of 
Civil  Defense,  Department  of  the  Army  under  Contract  DAHC  20-70-C-0287 
with  the  Institute  for  Defense  Analyses. 


Edward  S.  Pearsall 
May  1971 
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A  LAGRANGE  MULTIPLIER  METHOD  FOR 
CERTAIN  CONSTRAINED  MIN-MAX  PROBLEMS 


I 

INTRODUCTION 

Constrained  min-max  problems  are  constant-sum  two-person  games  in 
which  the  maximising  player  enjoys  the  advantage  of  ..oving  last 
and  both  players  select  strategies  subject  to  separate  side  condi¬ 
tions.  Problems  of  this  general  form  are  natural  archetypes  for  a 
variety  of  defense  resource  allocation  models.  Probably  the  most 
common  defense  application  is  to  the  allocation  among  cities  of 
shelters  and  antiballistic  missiles.  The  objective  of  the  defense 
is  to  minimise  the  casualties  that  may  be  inflicted  by  an  enemy 
with  a  given  arsenal  of  ballistic  missiles  and  a  full  view  of  the 
defense  allocation.  Among  the  authors  who  have  formulated  such 
models  as  constrained  min-max  problems  are  Robert  Kupporman, 

Smith  r  10]  ,  Eisen  and  'Owen  [  7]  .  An  excellent  mathematical 
survey  of  min-max  problems  in  general  has  beer,  provided  by 
Danskin  [3J. 

Since  the  minimizing  player  is  deprived  of  the  means  to  imple¬ 
ment  a  mixed  strategy,  constrained  min-max  problems  do  not  always 
possess  equilibrium  pairs,  i.e.,  strategies  that  are  simultaneously 
optimal  against  each  other.  This  inconvenient  characteristic 
separates  min-max  problems  from  the  well-developed  body  of  theory 
for  games  with  sacdle-point  solutions.  Methods  designed  to  locate 
saddle-points  will  succeed  only  fortuitously  in  solving  min-mex 
problems . 

Equally  important,  the  absence  of  saddle-point  solutions  pre¬ 
vents  the  straightforward  use  of  Lagrange  multipliers  in  dealing 
with  the  players'  constraints ,  Nevertheless,  several  procedures 
employing  Lagrange  (multipliers  (or,  equivalently,  a  convex  duality 
theorem)  have  been  proposed  for  min-sax  problesas  with  constraints, 
a: a-  (2)  and  Danskin  [3]  prove  Lagrange  multiplier  theorems  in 


1 


connection  with  their  use  of  directional  derivatives.  However,  the 
constraints  pertain  only  to  the  minimizing  player  and  directional 
derivatives  typically  lead  to  cumbersome  processes  for  calculating 
solutions  in  applied  work.  A  more  workable  approach  is  the  "double" 
Lagrange  multiplier  method  proposed  by  Everett  [5]  and  Pugh  [9 3* 
Unfortunately,  the  circumstances  under  which  the  method  might 
yield  optimal  solutions  are  not  described  by  Everett  and  Pugh. 
Moreover,  their  method  is  clearly  capable  of  yielding  strategies 
that  are  not  optimal  solutions  to  constrained  min-max  problems. 

Pugh  attempts  to  deal  with  this  difficulty  by  offering  "verification" 
techniques  that,  in  effect,  provide  bounds  for  solutions.  However, 
these  bounds  are  of  doubtful  value  in  applied  work.  Most  recently, 
Owen  [7J  nas  employed  a  convex  duality  theorem  to  transform  a  con¬ 
strained  min-max  problem  into  a  pure  minimization  problem.  Unfor¬ 
tunately,  Owen  fails  to  exhibit  solution  strategies  for  the 
maximizing  player. 

In  this  paper,  we  shall  employ  Owen’s  procedure  to  obtain  a 
specialization  of  Everett  and  Pugh’s  double  Lagrange  multiplier 
method.  The  method  applies  to  problems  in  which  the  maximizing 
player’s  strategies  may  be  segregated  by  tactics  and  probabilistically 
mixed,  For  example,  in  attacking  a  city  defended  by  antiballistic 
missiles,  the  maximizing  player  might  employ  strategies  designed  to 
leak  through,  suppress  (by  destroying  the  radars),  or  exhaust  the 
defense— or  he  may  choose  not  to  attack  the  city  at  ail.  Normally, 
the  maximizing  player  would  r.ot  be  indifferent  between  two 
strategies  associated  with  the  same  tactic,  i.e.,  having  decided 
to,  say,  exhaust  the  defense,  the  maximizing  player  will  typically 
find  that  there  is  a  single  best  number  of  warheads  for  this  pur¬ 
pose.  On  the  other  hand,  a  clever  defense  may  leave  the  maximizing 
player  indifferent  between  the  best  strategies  associated  with 
two  or  more  distinct  tactics.  For  example,  antiballistic  missiles 
may  have  been  installed  in  the  city  in  numbers  just  sufficient  tc 
leave  the  maximizing  player  indifferent  between  exhausting  the 
defense  and  not  attacking  the  dry  at  all.  We  shall  allow  the 


maximizing  player  to  resolve  such  dilemmas  by  probabilistically 
mixing  strategies  associated  with  different  tactics. 

This  is  a  novel  formulation  for  constrained  min-max  problems 
and  requires  some  elaboration  which  we  shall  provide  in  Section  II. 
In  Section  III  a  convex  duality  theorem  is  employed  to  derive  a  set 
of  sufficient  conditions  for  solutions  to  a  class  of  constrained 
min-max  problems  with  mixed  strategies.  An  even  stronger  set  of 
sufficient  conditions  may  be  derived  if  certain  functions  are 
differentiable  and  quasi-convex.  These  conditions  are  derived  in 
Section  IV  and  provide  the  basis  for  the  presentation  of  the 
modified  double  Lagrange  multiplier  method  (Section  V).  In  the 
concluding  Sections  the  characteristics  of  constrained  min-max 
problems  to  which  the  method  is  best  suited  are  described  and  a 
simple  ABM/shelter  deployment  problem  is  solved  to  provide  an 
illustration . 


II 


THE  CONSTRAINED  MIN-MAX  PROBLEM  (WITH  MIXED  STRATEGIES) 

By  convention,  the  outcome  of  a  constant -sum  two -person  game  is 
measured  as  a  payoff  to  one  of  the  players  from  which  the  other  ’s 
loss  differs,  at  most,  b,  a  constant  amount.  If  we  let  X  and  Y 
denote  strategies  for  the  maximizing  and  minimizing  players,  re¬ 
spectively,  then  this  payoff  (to  the  maximizing  player)  may  be 
represented  by  the  real-valued  function  H(X,Y).  The  constraints 
on  the  players '  uses  of  resources  may  be  represented  by  the  inequality 
systems  C(X)  <  c  and  B(Y)  1  b  where  C(X)  and  B(Y)  denote  real-valued 
vector  functions  and  c  and  b  are  vectors  of  the  available  quantities 
of  the  players’  resources.  Moreover,  the  players  are  restricted  in 
their  selection  of  strategies  to  the  elements  X  and  Y  of  the  sets 
S  and  T.  The  function  H  is  defined  over  the  product  set  S  x  T,  the 
vector  function  C  is  defined  over  S,  and  B  is  defined  over  T. 

(Throughout  the  paper  vectors  are  taken  to  be  row  vectors  unless 
marked  with  an  apostrophe  to  indicate  the  transpose.) 

A  general  statement  of  the  constrained  min-max  problem  is  to 
find  a  strategy  X°  for  the  maximizing  player  and  a  strategy  Y° 
for  the  minimizing  player  that  satisfy  the  requirements: 

1.  The  strategies  are  feasible  for  boch  players:  X°€S,  C(X°)  1  c 
and  Y°€T,  B(Y°)  <.  b. 

2.  The  strategy  X°  maximizes  the  payoff  function  given  the  strategy 
Y°  subject  to  X€S  and  C(X)  <  c: 

H(X°,Y°)  >  H(X,Y°)  for  all  X€S,  C(X)  <  c. 

3.  The  strategy  Y°  minimises  the  maximum  payoff  in  X€S,  C(X)  ^  c 
subject  to  Y€T  and  B(Y)  <  b: 

|  H(X,Y)  |  X€S,  C(X) 


H(X°,Y°)  <  Sufy 


*  C \ 


for  all  Y€T,  B(Y)  <  b. 


A 


(Note  that  a  constrained  max-min  problem  may  always  be  converted  to 
a  min-max  problem  with  the  same  constraints  simply  by  reversing  the 
sign  of  the  payoff  function.) 

Mixed  strategies  are  not  conventionally  believed  to  serve  any 
useful  purpose  in  the  mathematical  theory  of  min-max  problems.  The 
minimicing  player  cannot  implement  probabilistic  mixes  of  strategies 
and  the  maximizing  player  appears  to  gain  nothing  by  using  them. 

In  short,  it  is  impossible  to  visualize  any  useful  function  of  mixed 
strategies  in  an  actual  play  of  the  game.  From  this  observation  it 
is  easy  to  presume  that  mixed  strategies  for  the  maximizing  player 
can  have  no  role  in  the  mathematical  treatment  of  min-max  problems. 
This  has  been  the  presumption  of  virtually  all  previous  work.  (Note, 
however,  that  Pugh  [B]  and  Owen  [7]  in  their  work  replace  the  payoff 
function  with  its  concave  envelope.  This  device  implicitly  admits 
mixed  strategies  since  there  may  be  points  on  the  concave  envelope 
that  correspond  only  to  weighted  averages  of  points  on  the  real 
payoff  function.  These  payoffs  can  be  achieved  only  by  probabilis¬ 
tically  mixing  the  maximizing  player’s  strategies.) 

In  standard  game  theory  the  introduction  of  mixed  strategies  for 
both  players  serves  to  insure  the  existence  of  a  saddle-point  solu¬ 
tion.  That  is,  there  exist  mixed  strategies  that  constitute  an 
equilibrium  pair.  The  introduction  of  mixed  strategies  for  just 
the  maximising  player  is  not  alone  sufficient.  The  resultant  game 
will  not  possess  a  saddle-point  unless  the  min-max  problem  had  an 
equilibrium  pair  of  strategies  in  the  first  place. 

Nevertheless,  mixed  strategies  for  che  maximizing  player  may  be 
made  to  serve  a  useful  purpose  by  the  interesting  expedient  of  re¬ 
defining  the  concept  of  an  equilibrium  pair  and  then  exhibiting 
a  correspondence  between  such  redefined  equilibrium  pairs  and 
solutions  to  constrained  min-max  problems.  A  slightly  modified 
version  of  the  ’’seesaw"  problem  provides  a  simple  illustration  of 
this  use  of  mixed  strategies.  The  maximizing  player  chooses  a  point 
between  the  right  end,  X=l,  and  the  left  end  X=-l,  of  a  seesaw  but 
he  may  not  sit  directly  over  the  fulcrum,  i.e.,  -  1  <  X  1  1  and  X  +  0. 
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The  minimizing  player  chooses  the  angle  of  the  seesaw,  -  n/2  <  Y  <  n/2. 
The  payoff  is  the  height  of  the  maximizing  player’s  seat  above  the 
fulcrum,  H(X,Y)  =  X  sin  Y.  This  problem  has  no  equilibrium  pair  of 
strategies  (X°,Y°).  However,  a  mixed  strategy  may  be  used  to 
demonstrate  that  V°  =  0  is  optimal  for  the  minimizing  player.  Two 
tactics  may  be  descerned  for  the  maximizing  player  by  observing  that 
either  X=-lorX=lis  optimal  for  the  maximizing  player  for  any 
Y  (both  are  optimal  if  Y  =  0).  Let  H£(Y)  ana  H^(Y)  be  the  payoffs 
that  result  from  a  choice  of  the  left  and  right  ends  of  the  seesaw, 
respectively.  Then,  H*(Y)  =  -  sin  Y  and  H*(Y)  =  sin  Y.  A  probabil¬ 
istic  mix  of  these  two  strategies  yields  the  expected  payoff  u(-sin  Y) 

+  (1.  -  u)  sin  Y  =  (1  -  2u)  sin  Y  where  0  <  u  <  1.  The  probability 
u°  =  .5  and  the  strategy  Y°  =  0  constitute  an  equilibrium  pair  in 
the  restricted  sense  that : 

(1  -  2u°)  sin  v°  :>  (i  -  2u)  sin  Y°  for  all  0  S  u  <  1 

and 

(1  -  2u°)  sin  Y°  <  (1  -  2u°)  sin  Y  for  all  -  H/2  ^  Y  %  n/2. 

The  maximizing  player’s  mixed  strategy  is  to  sit  on  the  right  end  of 
the  seesaw  (X°  =  1)  with  probability .5  (u°  =  .5)  and  to  sit  on  the 
opposite  end  of  the  seesaw  (X°  =  -1)  with  the  same  probability  .5 
(1  -  u°  =  .5).  That  =  0  satisfies  Statement  3  of  the  definition 
of  a  solution  to  the  min-max  problem  follows  from  the  fact  that : 

(1  -  2u°)  sin  Y  <  Kaxx  jx  sin  Y  1 -1  i  X  1  1,  X  /  o} 

for  all  -  n/2  %  Y  %  n/2. 

Either  X°  =  -  1  or  X°  =  1  satisfy  Statement  2  and  both  (X°  =  -  1, 
y°  =  0)  and  (X°  =1,  Y°  =0)  satisfy  the  feasibility  requirement. 
Statement  1.  Now,  the  mixed  strategy  we  have  derived  for  the 
maximizing  player  is  completely  superfluous  so  far  as  playing  the 
seesaw  game  is  concerned.  Its  function  in  the  argument  above  is 
solely  to  permit  us  to  identify  an  optimal  strategy  for  the 
minimizing  player.  This  approach  to  min-max  problems  will  be 


6 


formalized  and  extended  to  problems  with  constraints  in  the  remainder 
of  the  paper  (The  solution  to  the  seesaw  problem  is  derived  with 
directional  derivatives  in  Danskin  [3].) 

For  the  special  class  of  problems  dealt  with  in  this  paper,  we 
assume  that  the  maximizing  player’s  strategies  may  be  assembled  by 
tactics  and  that  strategies  associated  with  different  tactics  may 
be  probabilistically  mixed.  Let  j  =  l,...,n  be  an  index  of  tactics. 
For  each  tactic  j  we  denote: 


sj 

w 

uj 


a  set  of  strategies  Xj 
a  payoff  function 

a  resource  vector  function  for  the  maximizing  player, 
a  probability  of  use. 


We  also  assume: 

1.  The  strategies  are  representable  as  nonnegative  row 

vectors  of  real  numbers,  i.e.,  =  {  Xj  ^  ^  I  Xj  %  °}  • 

2.  The  resource  vector  functions  are  linear,  i.e.,  Cj(Xj)  =  X^A^ 

where  A.  is  a  real  matrix  of  appropriate  dimensions. 

3.  The  systems  X^  >  0,  X^  l  c  all  have  at  least  one  solution. 

4.  The  functions  H.t(Xj,Y)  are  finite-valued  and  concave  in  X^  >  0 
for  all  Y€T*  B(Y)  %  b. 

5.  The  functions  H^(V,Y)  defined  by: 

H|(V,Y)  =  Supx  {  Hj(X^,Y)  -  Vn'X'  |  X.  >  o} 

are  continuous  in  V  >  0  [A^  and  denote  the  transposes  of 
k.  and  Xy  ]  for  all  Y€T,  B(Y)  %  b. 

6 .  c  >  0 . 

The  purpose  in  assembling  strategies  by  tactics  is  to  provide 
functions  H|(V,Y)  j=l,...,n  with  convenient  mathematical  properties 
when  Y  is  defined  as  a  nonnegative  real  vector.  Specifically,  it  is 
generally  necessary  to  define  tactics  in  such  a  way  that  the  functions 
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Hj(V,Y)  j  l,...,n  are  continuous  and  differentiable  for  V  >  0  and 
V  >  0.  Sorting  the  maximizing  player’s  strategies  by  tactics  usually 
does  little  conceptual  violence  to  a  constrained  min-max  problem. 
Indeed,  useable  divisions  of  the  strategy  set  S  are  often  suggested  by 
the  structure  of  the  problem  itself  and  probabilistically  mixing  the 
strategies  associated  with  different  tactics  causes  few  difficulties 
in  the  interpretation  of  solutions. 

The  constrained  min-max  problem  with  mixed  strategies  is  to  find 
strate9ies  X.  j=l» . . . »n  for  each  of  the  maximizing  player’s  tactics, 
probabilities  u^.  3=1,..., n,  and  a  strategy  Y°  for  the  minimizing 
player  such  that : 


1. 


2. 


X°  ^  0  j=l,..., n,  u°  >  0  j=l, 
and  Y°€T,  B(Y°)  1  b. 


n 

„  o 
3=1  ^ 


=  1, 


n 

E 

3=1 


uSxSAj  =  0 


"  uV(X°  ■?)  >  Z  u.H.CX-.Y0)  for  all  X.  >  0 
j=l  3  3  3  -  j=1  3  3  3  3  = 


3=1, 


,n, 


i. 


n 

u.  ^  0  j=l,..., n,  T.  u.  =  1,  ?_  u.X.A.  <  c. 

3=1  3  j=l  J  2  3  ~ 


Jj.  S  s“PxjlU.  {  |  X.  >  0 


n  n  \ 

u.  ^  0  j=l,..., n,  I  u.  =  1,  I  u.X.A.  <  c  v 
3  j=i  3  3=1  3  3  3  "  f 

for  all  YCT,  B(Y)  <  b. 

xhis  definition  is  an  obvious  analogue  of  the  definition  of  a  solu¬ 
tion  provided  earlier  for  constrained  min-max  problems  without  mixed 
strategies . 

The  definition  incorporates  expectational  statements  of  the  payoff 
n 

function,  7  u-H.  (X.  ,Y),  and  the  maximizing  player’s  resource 
3=1  333 

n 

constraints,  Z  u^X.A.  <  c.  That  is,  the  solution  to  the  constrained 
j-l  333- 

min-max  problem  with  mixed  strategies  establishes  a  random  process  in 
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which  the  strategies  X°  j=l,...,n  are  played  with  probabilities 
Uj  j=l,...,n.  The  outcome  of  the  game  is  the  expected  payoff 

T.  p°H.(X°,Y°)  and  the  expected  use  of  resources  by  the  maximizing 
j=l  3  3  3 
n 

player  E  u  .X°A .  does  not  exceed  the  vector  of  available  quantities 
j=l  3  1  3 

c.  However,  there  is  no  guarantee  that  the  expected  payoff  will  be 
achieved  and  that  the  maximizing  player’s  resource  constraints  will 


be  respected  in  any  single  play  of  the  game.  In  most  applications, 


therefore,  it  would  be  preferable  to  solve  for  strategies  X^  j-1,. 

such  that  if  >  0  then  X^A.  1  c  and  H.(X^,Y°)  -  E  p^H-CX^jY0). 

J  J  J  J  J  j  —  2,  J  J  * 


With  this  condition  observed  any  play  of  the  game  must  result  in  the 
same  actual  payoff  at  a  cost  in  resources  to  the  maximizing  player 
that  never  exceeds  the  amounts  available.  Unfortunately,  the  in¬ 
clusion  of  this  additional  restriction  in  the  definition  of  a  solution 
makes  mixed  strategics  superfluous  and  leaves  us  with  a  mathematically 


less  tractable  problem. 

A  sensible  second-best  procedure,  then,  is  to  solve  for  strategies 
that  satisfy  only  Statements  1  through  3  of  our  definition.  Whether 
or  not  this  procedure  produces  acceptable  approximations  can  be 
determined  after  the  fact .  Variations  in  payoff  and  resource  use  may 
be  computed  and  if  they  are  small  relative  to  the  expected  values, 
little  distortion  is  likely  to  result  from  accepting  them  as  approxi- 

mations.  In  any  case  the  expected  payoff,  ^  provides  a 

useful  upper  bound  to  the  solution  value  of  the  payoff  function  if  we 


add  the  condition  X°A.  <  c  and  K.(X^,Y°) 

j  j  ~  j  j 

for  all  tactics  j--l,...,n. 


7.  Y°)  if  u*  >  C 

j=l  1  J 


Q 


Ill 


SUFFICIENT  CONDITIONS  FOR  A  SOLUTION 

Lee  V  denote  a  vector  of  nonnegative  real  numbers  whose  elements 
correspond  to  the  elements  of  the  resource  vector  function 
n 

r.  Sufficient  conditions  for  a  solution  to  the  constrained 

j=l  3  3  3 

min-max  problem  with  mixed  strategies  are: 

1.  X°.  >  c  j=l, . . .  ,n,  u?  >  0  j=l,...,n,  ?  U?  =  1,  Y°€T,  V°  >  0. 


2a.  H ^ ( X j )  -V°A'X°'  ^  H^X^Y0)  -V°A^X'  for  all  Xj  ^  0  j=l, 

2o.  ^  a?  [H.(X^,Y°)  -V°A'X^]  >  V  u.[H.(X°  Y°)  -  vVx^'] 

j=a  J  3  j  J  ~  j=i  j  j  3 


for  all  u-  >  0  j=l,...,n,  7.  u.  =  1. 


n  0  n  T  ^ 

2c.  V°|  “  u°X°A.  -  c  >  V|  v  U°X°A.  -  c  for  all  V  >  0. 

lj=l  3  3  3  J  Lj-1  3  3  3  J 

n  r* 

3a.  T  u°H* (V° , Y° )  +  V°c'  <  7  u°H:;(V,Y)  +  Vc' 


1=1  3  3 


=  3=1  3  3 

for  all  V  >  0,  Y€T,  3(Y)  <  b. 


3b.  3(Y  )  <  b 


where  the  functions  Hv  are  defined  as: 


Hj(V,Y)  =  Supx  |  Hj(Xj ,  Y)  -  VA'X^Xj  >  oj 

and  C  j'  denotes  the  transpose  of  a  matrix.  The  numbering  of  the 
conditions  has  been  chosen  to  roughly  correspond  to  the  numbered 
statements  in  the  definition  of  a  solution. 
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Condition  1  largely  repeats  Statement  1  of  the  definition.  Non¬ 
negative  conditions  on  the  solution  value  of  the  vector  V  have  been 

n 

added  and  the  inequalities  7  u°X°A .  <  c  and  B(Y°)  <  b  have  been 

j=l  :  3  3  " 

deleted.  The  reasons  for  these  alterations  will  become  apparent  as 
we  proceed  to  show  the  origins  of  Conditions  2  and  3. 

Condition  2,  together  with  the  relevant  components  of  Condition  1 
.nay  be  used  to  show  that  Statement  2  of  the  definition  of  a  solution 

n 

and  the  inequality  T.  u°X°A.  <  c  are  satisfied.  For  this  demon- 

j=I  3  3  3 

stration  we  require  the  Saddle-Point  Theorem  of  mathematical  program¬ 
ming  .  For  ovr  purposes  a  convenient  statement  of  the  Theorem  is 
the  following:  Let  F(X)  be  a  real-valued  function  and  G(X)  a  real¬ 
valued  vector  function  of  the  strategies  X  in  S.  If  the  strategy 
X°€S  and  a  real  vector  V°  i  0  constitute  a  saddle-point  of  the 
Lagrangian  function  F(X)  -  V[C(X)]'  xn  X€S  and  V  £  0»  i.e., 

F(X°>  -  V°rG(X0)]'  >  F(X)  -  V°[G(X)]'  for  all  X6S 
and 

F(X°)  -  V°[G(X°)r  <  F(X°)  -  V[G(X°)J#  for  all  V  *  0 
then  X°  maximizes  i ' ( X )  in  S  subject  to  G(X)  <  0,  i.e.» 

F(X°)  >  F(X)  for  all  X€S,  G(X)  <  0 

and  G(X°)  <  0.  The  Theorem  can  be  found  in  most  texts  on  mathematical 
programming  such  as  Karlin  [6].  The  proof  is  not  dependent  on  the 
characteristics  of  the  set  S. 

To  apply  the  Theorem  let : 

S  =  |  X^,u^  j=l» . .  •  »n  j  Xj  \  0  j=l, . . .  ,n,  Uj  >0  jsi» . . .  »n, 
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n 


F(X)  =  7  d,H.(X.,Y°) 
3=1  3  3  3 


and 


n 


G(X)  =  2  u.X.A.  -  c. 
j=l  3  3  3 

Condicions  2a  and  2b  may  be  used  to  construct: 

2  U°H.(X°,Y°)  -  v°[  2  u°X°A.  -  c y 


5=1 


3  3'  3 ' 


=i  ]  ]  3 


j=l 


n  n 

>  2  d -H  .(X  . , Y  )  -  V°[  7  u,X,A,  -  c]' 

~  j=l  3  3  3 


j=l  3  3  ^ 


for  all  X.  >  0  j=l,..., n,  u.  >  0  j=l,...,n,  v  u  -  i 
3  -  J 


From  Condition  2c : 


2  u°H  (X°Y°)  -  V°[  2  u°X°A.  -  c]' 
5=1  3  3  3  j=l  3  3  3 

-  d°H  (X°,Y°)  -  V[  " 

3=1  3  3  -  3  = 


n  n 

<  T.  ^H.(X°,Y°)  -  V[  2  u?X?A_-  -  c  ]'  for  all  V  >  0. 
-  j=1  J]  ]  3  =!  3  3  J 


Since 


n 


X?  >  0  j=l,  —  ,r. ,  >  0  j=l,..., r.,  r  u°  =  1  and  V°  >  0 

J  J  ~  ~  — 1  J  ~ 

j  —X 

we  must  have  a  saddle-point.  By  the  Saddle-Pcir.t  Theorem: 
n  n 

I  u^H.(X^,Y°)  >  ^  u.H.(X_.,Y°)  for  all  X.  >  0  j=l,...,n, 

3=1  J  3  3  j=l  3  3  J  3 

n  n 

u  •  >  0  j=l, . . .  ,n,  7  u  .  =  1,  2  u.X.A.  <  c 
3  "  3=1  3  5=1  3  3  3  ~ 

which  is  Statement  2  of  the  definition  of  a  solution  and 


n 

v  u.X°A.  <  c. 
j=l  3  3  3  ~ 
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Xl. 


The  first  step  of  the  demonstration  that  Statement  3  of  the 
definition  of  a  solution  is  implied  by  Conditions  1  through  3  is  to 


n 

exhibit  a  saddle-point  of  the  quantity  Z  u.tf$(V,Y)  +  VC#. 

j=l  3  3 


Conditions  3a  and  3b: 


From 


n  n 

v  U°H?(V°,Y°)  +  V°c'  <  r  a°H*(V,Y)  +  Vc' 
j=i  3  J  ~  3=1  3  3 

for  all  YcT,  B(Y)  %  b,  V  >  0 

and  B^)  ^  b. 

Condition  2b  provides: 

”  u°H*(V°,Y°)  +  V°c'  >  r  u.H*(V°,Y°)  +  V°c' 

j=l  3  3  "  j=l  3  3 


n 

for  all  u.  >  0  j=l,...,n,  T  u •  =  1. 

3  ~  j=l  3 

n 

So  we  apparently  have  a  saddle-point  of  the  expression  T  p  .H^(V,Y) 

j=l  33 

+  Vc'  where  the  minimization  proceeds  over  Y€T,  B(Y)  <  b,  V  >  0  and 

n 

the  maximisation  is  taken  over  u.  >  0  j=i . .  X  u4  =  1.  This  is 

3  "  j=l  3 

the  equilibrium  pair  in  a  restricted  sense  that  we  referred  to  earlier 
Since  the  order  in  which  the  minimisation  and  maximization  operations 
are  performed  is  i material  w hen.  dealing  with  saddle-points  we  have: 

Z  u VvVC,VP)  ♦  vV  <  Hih.Wa*  \  r  U.H^CV.Y)  +  Vc' 
j=i  3  •’  “  Uj  3  3 


jv>  C,  ^  >  0  j=l,. . ,  ,n,  I  u-  =  it  for  all  YST,  B(Y)  <  b. 
3  "  j=i  3  ) 

The  next  3tep  is  to  show  that  Kir^Hax^  equals  Max^  Min^. 

This  is  done  by  appealing  to  the  familiar  Min-Max  Theorem  of 


game  theory.  The  following  statement  of  the  Theorem  approximates  that 
given  in  Karlin  [6J.  Let  F(X,Y)  be  a  real-valued  function  of  X  in  S  and 
Y  in  T  where  both  S  and  T  are  closed,  bounded  and  convex  sets  in  £n. 

If  F  is  continuous,  convex  in  Y  for  each  X,  and  concave  in  X  for  each 
Y,  then: 

MinyMaxx  |f(X,Y)|X€S,  Y€t|  =  MaxxMiny  |f(X,Y)IX€S,  Y6T  }  . 

To  apply  the  Min-Max  Theorem  define: 

S  =  <  o  .  j-1, . . .  ,n  j  u  .  >  0  j=l, . . .  ,n,  T  u  .  =  1  >  , 

I  3  I  3  -  j  ? 

T  =  j  V  |  V  ^  0,  Vc'  <  Max..  Jh^(0,Y)  -  H^(0,Y)j|, 
and 

n 

F(X,Y)  =  T.  U,H*(V,Y)  +  Vc\ 
j=l  3  J 

It  is  obvious  that  S  is  closed,  bounded  and  convex  and  that  T  is 
closed  and  convex  for  any  Y€T,  B(Y)  <  b.  Tc  show  that  T  is  also 
bounded  we  must  recall  that  c  >  0  and  that  the  functions 

Hj(Xj,Y)  j=l,...»n  are  finite  valued  in  X^  >  0  j=l n  for  Y€T» 

B(Y)  <  b.  Then  Max.  1  H*(G,Y)  -  H_.(Q,Y)\  <  »  and  the  set  T  must  be 

=  3  13  .»  * 


bounded.  The  function  u-H*(V,Y)  *  Vc*  is  clearly  continuous  and 

3=1  3  3 

concave  in  the  variables  j*l,,..,:t  for  all  values  of  the  vector  V. 
The  function  is  continuous  in  V  tor  any  u.-  j=l»...»n  because  the 

J 

functions  H?(V,y)  j=l,...,n  are  each  assumed  to  be  continuous  in  V. 

3 


Since-  Uj  <► 


3=1.  . 


•  hC  • 


j  d*?Y,Y)  ♦  Vc'  is  convex  in 
jj 


V  if  the  functions  H*(V,V)  are  each  convex  functions  of  V. 

Choose  Vx  £  0  and  v/>  0,  the:,  for  an  0  <  t  <  1: 


fIjCcVl  +  Cl-t)V2>Y)  =  Supx  | H  ^ ( Xj , Y)  -  (tV1+(l-t)V2)A'X'|Xj  :>  0 

i  SupXj{tHj(Xj»Y)  -  cVlAjXj|Xj  ^  °}  +  SuPx. 

-  U-tOVjA'X^X..  >  o}  =  tH*(V1,Y)  +  (1-t)  HJ(V2,Y) 
which  proves  convexity. 

Since  the  conditions  of  the  Min-Max  Theorem  are  met  for  all 
YET ,  B(Y)  g  b:  : 

(  n  .  I  n 

Minv  Max  {  r  u^(V,Y)  +  Vc'  U.  >  0  j=l . .  V  w  =  1, 

j  lj=l  J  J  D  j  j=l  3 

V  i  0,  Vc'  g  Max..jH*C0,Y)  -  H^O.YjU 

in  n 

=  Ma* .  Minv  [  Z  U-H*(V,Y)  +  Vc'ja..  >  0  j=l,...,n,  I  u .  =  1, 

Jj  (j=l  3  J  1  -  "  v  j=l  3 

V  ^  0,  Vc'  1  MaXj|Hv(0, Y)  -  H-j(°»Y)U. 

The  bound  Vc'  g  Max^.|Hv(0,Y)  -  H^(0,Y)|  is  superfluous.  Consider 

any  saddle-pointCu^  j=l,...,n,  V).  At  such  a  point  we  must  have: 

n  n  ; 

I  tI.K*(V,Y)  +  Vc'  <  Z  u -H?(V,Y)  +  Vc' 

3=1  3  3  ~  j=l  3  J 

for  any  V  ^  0,  Vc'  g  Max^HnCO/Y)  -  H.j(0,Y)}  . 

Since  the  inequality  must  hold  at  V  =  D: 

n  n 

v  J  Hv(V,Y)  +  Vc'  <  Z  G.Hv(0,Y). 
j=l  3  3  "  j=l  3  3 


A  L.l  O  y 


H* ( V, Y)  =  Supx  I^.CX^Y)  -  VA'X^|X.  :>  o}  ^  H..(0,Y)  j=l,...,n. 
Vie  may  substitute  and  rearrange  terms  to  obtain  the  inequality: 


Vc'  <  V.  a,[H*(0,Y)  -  H.(0,Y)]. 
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n 

On  the  right-hand  side,  u  ■  ^  0  j=l,...,n  and  7  u..  =  1,  therefore, 
,  3=1'  3 

Vc  <  MaXjjHj (0, Y)  -  H^(0,Y)|  .  Thus  the  bound  Vc'  <  Max^.  |h*?(0,Y) 

Hj(0,Y)|  is  always  respected  and  its  inclusion  to  meet  the  condi¬ 
tions  of  the  Min-Max  Theorem  is  purely  formal.  Omitting  the  bound 
we  have: 


I  U?H*(V°,Y°)  +  vV  <  Me 


l  3  3 


Minv  7  H,H* 
1  (j=l  J  J 


(V, Y)  +  Vc' 


V  ^  °»  >  0  j=l,...,n,  I  u.  =  lj  for  all  Y€T,  B(Y)  <  b. 

J  3=1  J  !  - 

Our  next  step  is  to  show  that: 


Minv j  ^  PjH^CVjY)  +  Vc'|v  ^  o| 

=  Sup  {  7,  u-H.(X.,Y)|X.  >  0  j=l, 
J  ( j=l  3  3  3  1  J  " 


n 

7 

3=1 


f 

i 


n 

for  all  u .  >  0  j=l,...,n,  7  n .  =  1  and  YST,  B(Y)  <  b. 

3  ~  j=!  3 


Our  procedure 


Is  similar  to  that  followed  by  Owen  [7]  and  makes  use  of  the  same 


Convex  Duality  Theorem.  The  Theorem  is  due  to  Rockafellar  ["Si  who 


proves  it  in  a  somewhat  more  general  form  than  that  employed  by  Owen. 
Let  F(X)  be  a  concave  function  of  a  vector  X  in  En.  The  conjugate  of 
F  is  the  function  F*(U)  of  U  in  En  defined  by: 


F*(U)  =  Sup^FCX)  -  ax' |  . 

Let  G(Z )  be  a  convex  function  of  a  vector  Z  in  En.  The  conjugate  of 
G  is  the  function  G*(V)  of  V  in  En  defined  by: 

G*(V)  =  Infz|G(Z)  -  VZ'|  . 
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If  the  system  X  ^  0,  XA  ^  Z  (ft  is  3  real  matrix  of  appropriate 
dimensions)  has  a  solution  in  the  domains  (of  finiteness)  of  F  and 
G,  then: 

SupX,z{F(X)  ~  G(Z)  I  X  >  0,  XA  •g  zj 
=  Minu>vjF*(U)  -  G*(V)  |  V  £  0,  VA'  £  tl}  . 

To  apply  Rockaf ellar ’s  Theorem  we  define: 


and 


F(X) 


n 

S  UHH-(X,5Y)  if  X.  >  0  j=l . . 

j  — 2_  J  J  J  J  ~ 


-  ®  otherwise 


G(Z  )  = 


0 


CO 


if  Z  =  c 
otherwise . 


To  verify  that  the  conditions  of  the  theorem  are  met  we  observe, 
first,  that  F(X)  is  concave  in  X^  j=l,...,n  for  >  0  j=l,...,n 
provided  that  the  functions  H^(X^,Y)  are  concave  in  X^  ^  0  for  all 
Y6T,  B(Y)  1  b.  But  this  has  already  been  assumed.  Second,  the 
function  G(Z)  is  trivially  convex  in  Z.  Lastly  we  must  show  that 

n 

the  system  X.  >  0  j=l,...,n,  E  (i.X.A.  <  Z  has  a  solution  in  the 
3  j  =1  3  3  3 

set  of  points  for  which  F  and  G  are  finite  for  any  lu  :>  0  j=l,...,n, 
r. 

E  (j  .  =  1.  Since  G(Z)  is  finite  only  at  Z  =  c  this  condition  is 

j=l  3 


met  only  if  X.  ^  0,  XI.  l  c  is  a  solvable  system  and  H.(X_ 


.Y') 


finite -valued  for  at  least  one  solution  for  each  tactic  j=l,...,n. 
We  have  assumed  that  Xj  :>  0,  X^Aj  i  c  is  solvable  and  that  H^(Xj,Y) 


is  finite-valued  in  X^  >  0  for  all  \€T,  B(Y)  <  b  for  every  tactic 
j=l,...,n.  The  conjugate  functions  F*(U)  and  G*(V)  are: 
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fn  n 

E  UhH.(X.,Y)  -  £  U.X^I 

j=l  3  3  3  j=l  3  31 


Xj  ^  0  j  1 9  •  •  •  » 


n! 


and 


G*(V)  =  -  Vc'. 

By  the  Convex  Duality  Theorem: 
Sup 


(  n  n  | 

T  u.H.(X.,Y)  -  G(Z)|X.  >  0  j=l,...,n,  £  u-X.A,  <  Z, 

j,Z  (j=l  3  3  3  '  3  "  j=l  3  3  3  "  ! 

=  IUnuj,v  {s“Pxd  {  J,  Ujx^|xj  i  o  0=1,. ...nj 

+  Vc'|v  :>  0,  Vu-A^  |  U.  j=l,...,nl  for  all  p.  ^  0  j=l,...,n, 

J  J  J  /  J 


n 


£  u .  =  1  and  Y€T,  B(Y)  1  b. 

3=1  3 

Each  side  of  this  equation  can  be  simplified.  On  the  left-hand 
side  a  supremum  always  is  attained  at  Z  =  c.  On  the  right-hand  side 
the  supremum  within  the  brackets  is  a  nonincreasing  function  of  the 
vectors  U.  j=l,...,n.  Therefore,  a  minimum  can  always  be  found 
among  the  vectorr  V,  IK  j=l,...,n  such  that  Vu^A^  =  IK  j=l,...,n. 

On  the  right-hand  side  of  the  equation  we  may  omit  the  inequalities 

Vu  -h'.  >  U.  i=l,...,n  and  make  the  substitution: 

33=  3 

/  n  n  ) 

SupY  \  £  |i-H.(X.,Y)  -  £  U.XtlX.  >  0  j=l,...,n> 

X3  (3=1  3  3  3  3=1  3  3  3  » 

(  n  n  ) 

=  SupY  l  £  q-H.(X.,Y)  -  £  q.VA^X.,  IX,  ^  0  j=l, . . .  ,n  V  . 

Xj  1  3=1  3  3  3  3=1  3  3  3'  3  ) 

This  can  be  simplified  still  further  by  inserting  the  functions 
Hj(V,Y)  j=l, . . . ,n : 


18 


Su%  |  ji  lWVY>  '  X  i ° 


n 

=  7.  u,H;r(V,Y). 
j=l  3  3 


After  all  of  these  simplifications  have  been  incorporated  the 
equation  provided  by  the  Oonvex  Duality  Theorem  becomes: 

SupX.  j  .S.  UjHi(Xi,Y)|Xi  ^  0  3=l»...,n,  7  u-X-A.  <  c\ 

J  (  J=1  J  J  J  13  j=i  3  3  D  -  j 

V|  MjHJCV^)  +  Vc'|v  ^  oj 


=  Min 


for  all 


n 


a-?  ^  0  7  m  •  =  1  and  YcT,  B(Y)  <  b. 

J  j=l  J 

Combining  this  result  with  that  of  the  previous  step  yields: 


n 


A  '7HJCV°^)  +  vV  5  SupxjlUj  J X  WVY) 

•  n  n  ) 

|a  •  >  0  j-l,...,n,  u.  >  0  j=l,...,n,  I  (j .  =1,  £  u-X.A  <  c> 

J  "  j=l  3  j=l  3  3  3  =  J 

for  all  Y€T,  3(Y)  ^  b. 

iO  obtain  Statement  3  of  the  definition  of  a  solution  it  remains 
only  to  show  that : 


n 


T,  u^Hq(V°,YD)  +  V°c'  =  I  u°H.(X°,Y°). 
j=l  J  J  j=l  3  3  3 

Fro,:.  Condition  2 a  we  obtain: 

=  Hj(X?,y°)  -  V°ft'X°'  j=l . . 
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Multiplying  by  the  probabilities  j=l»...,n,  summing  and  adding 
o  / 

V  c  to  each  side: 


n 


E  u?H*(V°,YC)  +  V°c'  =  V  uuH  .(X^ ,Y°;  -  V°  [  r  u?X°A.  -  c] 


j=l  3  J 


3=1 


3  3'  3' 


3=1 


3  3  3 


From  V  >  0,  Condition  2c  and  the  previously  derived  inequality 
n 

I.  n -X°A .  <  c  we  may  infer: 
j=l  3  3  3  - 


Vu  [  ~  u°X°A.  -  cY  0. 
j=l  J  3  3 


Statement  3  of  the  definition  of  a  solution  to  the  constrained  min 
max  problem  follows  since  B(Y°)  <  b  is  simply  Condition  3b. 


fC 


IV 


THE  SUFFICIENT  CONDITIONS  WITH  DIFFERENTIABLE  QUASI -CONVEX 

FUNCTIONS 


A  somewhat  stronger  set  of  sufficient  conditions  may  be  derived  for 
problems  with  the  characteristics  (in  addition  to  those  already 
assumed — see  page  7. ) 


7. 

T  =  Iy  in  En  |  Y  >  o} 

8. 

H.(X.,Y)  j=l,..., n  are  differentiable  in  X. 

9  9  9 

^  0  j=l, . .  .,n 

and  Y  >  0 . 

9. 

Hj(V,Y)  3=1,..., n  are  differentiable  in  V  | 

0  and  Y  £  0. 

10. 

3(Y)  is  differentiable  and  quasi-convex  in 

Y  ^  0. 

If  Assumptions  ?  through  10  are  met,  the  sufficient  conditions  of 
the  previous  section  may  be  replaced  by  the  following  stronger  set. 
Let  W  denote  a  vector  of  nonnegative  real  numbers  whose  elements 
correspond  to  the  elements  of  the  resource  vector  function  B(Y) : 


1.  X~  :>  0  j=l, . 
W°  >  0. 


n 

..,n,  u°  >  0  j=l,...,r,  I  =  1,  Y°  >  0,  V°  ^  0 

J  "  3=1  J 


2a.  For  every  tactic  j  ar.c  eacu  element  of  : 


5H.(X^Y°) 

— TT7 - v  ajk  ^  0  anc 

Jk 

5H,(X°,Y°) 


"“mCTT 


-  V°a'  =0  if  X".  >  0, 


3k 


jK 


jK 


th 


where  n  denotes  the  k  row  of  A  . . 
jk  J 


?1 


2b.  For  every  tactic  j: 

£  U?CH.(X°  Y°)  -  vVx°']  ^  H.(X°  Y°)  -  vV'X?'  and 
j  =1  J  J  j  J  J  3  3  3  3 

r  u^CH.CX^,^)  -  vVx0']  =  H.CX°,Y°)  -  V°AtX°'  if  u? 
j-^JJJ  3  3  33  33  3 


n  n 

2c.  V°[  v  ut°x?A.  -  c]'  >  V  [  r  U°X°A .  -  c]'  for  all  V 

j=l  3  3  J  ~  j=i  2  3  3 


3a.  For  each  element  of  Y: 


n  0bH.<X°,X°) 

-  u<  — 1 w — 


3=1  11 


A  > 


+  w 


+  w 


SB(Y°) 


"59; 


>  0  and 


5B(Y°) 


“5y 


]■ 


if  YJ  >  0. 


3b.  W°[B(Y°)  -  b]'  >  W[B(Y°)  -  b]'  for  all  W  >  0. 


4.  Z  u°H-:(V,Y)  +  Vc'  is  quasi-convex  in  V  >  0  and  Y  >  0 
3=1  3  3 

one  of  che  following : 
a.  Z  u°X°A.  -  c  f  0 

4  — -j  J  3  3  :'- 

-i 


&H,(X°  Y°) 


*>•  .1  “j-k2 

3  =1  J  k 


>  0  for  at  least  one  element  Y, 


c . 


n  &H-(X;,Y°) 

F  u  -  — ’‘TKf - <  0  for  at  least  one  element  Y 


iY, 


j=l 

and  there  exists  a  Y  ^  C,  B(Y)  i  b  such  that  Y^  >  0. 
n  &H  (X°  Y°) 

d.  Z  —  qy -  ^  3  for  ac  least  one  element  Yt 


3=1  3 


and  the  functions  Hv(V,Y)  j=l,...,n  are  twice  differen 


in  the  neighborhood  of  V°,Y°. 


•  0. 

>  o. 


and 

of  Y. 

of  Y 

of  Y 
tiable 
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If  the  functions  H;j'(V,Y)  j=l,...,n  are  all  convex  in  V>  0  and  Y>  0, 
Condition  4  may  be  omitted. 

In  Condition  1,  Y°€T  has  been  replaced  by  Y°  >  0  and  the  non¬ 
negativity  condition  W°  >  0  has  been  added.  Conditions  2b  and  2c 
are  unchanged.  The  new  Condition  2a  may  be  readily  derived  by 
applying  a  familiar  lemma.  Let  F(X)  be  a  differentiable  concave 
function  of  a  vector  X  in  the  region  X  ^  0.  Let  X°  >  0  be  chosen  such 
that  for  each  element  Xk  of  X: 

^-0  “*£>«• 

Then  X°  maximizes  F(X)  in  X  >  0.  To  obtain  our  nev;  Condition  2a,  let: 


F(X)  =  H.CX.,^)  -  V°A'X'  . 

According  to  the  lemma  X°  maximizes  Hj(X^,Y°)  -  V°A^X^  >  0,  i.e., 

H.(X°,YC)  -  V°A'X°'  >  H.(X.,Y°)  -  V°A'X'  for  all  X .  >  0. 

Since  this  must  hold  for  all  tactics  j=l,...,n,  the  new  Condition  2a 
implies  the  old  one. 

To  derive  the  remaining  conditions  (and  the  inequality  W°  >  0) 
we  must  first  evaluate  the  derivatives  of  the  functions  H?(V,Y) 
j=l,...,n  in  the  region  of  (V°,Y°).  This  is  done  by  making  use  of 
an  inference  from  Condition  2a.  For  every  tactic  j  and  each  element 

x.k  of  Xj : 


«h.(x<a°) 


J 


Xjk  =  0* 


Take  the  total  derivative  of  H*(V,Y)  at  (V°,Y°): 
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&H  (X*  Y°)  j 

where  — *  denotes  a  vector  of  partial  derivatives  with 

1 

respect  to  the  elements  of  Y.  It  is  clear  that  the  maximum5 is 


bH.CX^Y0) 


preserved  when  dX.^  is  nonzero  only  if  — 

jk 

hH  .  (  ,  y®  ) 

Therefore,  dHv  =  — 2__i! - dY'-  A'X°'dV', 

hH^CV0,^)  0H.(X°  Y°) 


"  V°3jk  ■  j0" 


j  — 1 j • • • )H« 


bH*(V°,Y°) 

3  bv  XjAj  j=l>...,n. 

The  nonnegativity  condition  W°  ^  3  added  tc  Condition  1  and  our 
r.ew  Conditions  3  and  4  all  proceed  from  a  Theorem  by  Arrow  and 
Er.thoven  [1],  The  Theorem  is  stated  here  in  a  form  appropriate  for 
constrained  minima  and  quasi-convex  functions.  Let  F(X)  be  a 
differentiable  quasi-convex  function  of  the  vector  X,  and  let  G(X) 
be  a  dif terentiable  quasi-convex  vector  function,  both  defined  for 
X  0.  Let  X°  and  X°  (a  vector)  satisfy  the  Kuhn-Tucker  Conditions, 
specifically: 

X°  >  0,  \°  >  0 


0F(X 

oT 


sg(x  )! 

ox 


>°  [G(X°)]#  =  C, 


where  denotes  a  vector  and 


a  matrix  of  derivatives  with 


respect  to  the  elements  of  X,  and  let  one  of  the  following  conditions 
be  satisfied: 


Arrow  and  Enthoven’s  Theorem  let: 

n 

F(X)  =  7.  U?H^(V,Y)  +  Vc' 
j=l  3  3 

and  G(X)  =  3(Y)  -  b.  Both  F(X)  and  G(X)  are  assumed  to  be  quasi-convex 
in  v  >  0  and  Y  >  0 .  The  Kuhn-Tucker  Conditions  for  a  constrained 
minimum  are: 


V°  >  C,  Y°  ^  0,  W°  >  0 


WC  [3(Y°)  -  by  =  0. 

Substituting  for  the  derivatives  of  Hv(V,Y)  j=l, _ ,n  the  Kuhn-Tucker 

Conditions  may  be  rewritten  as: 

V°  ^  0,  Y°  :>  G,  W°  ^  0 

7  u°X°A .  -  c  <  0  and  V°  T  Z  u?X?A .  -  cl  =0 
3=1  333  Lj=1  3  3  3  J 
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n  „  M,(X°  V°) 


aiid 


u  — 


=  0 


W°CB(Y°)  -  b]'  =  0.  j 

Conditions  1  through  3  of  our  new  set  imply  that  the  Kuhn-Tucker 
Conditions  for  a  constrained  minimum  are  satisfied.  Condition  1 
includes  the  nonnegativity  conditions  V°  ^  0,  Y°  ^  0.  Condition  2c 
is  equivalent  to: 

I  U°X°A.  -  c  <  0  and  V°[  7  u^X°A.  -  c ]'  =0. 
j=l  3  3  3  =  j=1  3  3  3 

Condition  3a  may  be  used  to  derive: 


3=1 


N4'- 


and  W°[E(Y°)  -  b]'  =  0  follows  from  Condition  3b. 

The  remaining  stipulations  of  Arrow  and  Enthoven’s  Theorem  apply 
co  the  constrained  minimization  problem  as  follows: 

n  n  bK-CX^V0) 

a .  Z  u°X°a'  -  c.  <  0  or  Z  — - >  0 

j=1  j  3  h]  h  j=1  { 

for  at  least  one  element  of  V  or  Y^  of  Y  [a^  denotes 
the  hth  row  of  A^]; 
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>  0  or 


OyO  # 

l  u.-X  .a,  .  -  c, 
j=1  J  J  hj  h 


n  o  &HtCX°  V3) 

AU°~X — <0 


tor  at  least  one  element  V^  of  V  or  Y^  of  Y  where  there 
exists  a  Y  >  0  with  YR  >  0  such  that  B(Y)  <  b; 

n  _  .  n  bH.CX^.Y0) 

"  UjX.a'  -  ck  f  0  or  Z  v°.  -  -W  -  ?  0 

j=l  3  3  3n  K  j=l  3  &Yk 

for  at  least  one  element  V^  of  V  or  Y^  of  Y  and  the 

functions  Hr?(V,Y)  j=l,...,n  are  twice  differentiable  in  the 

neighborhood  of  (V°,Y°); 


d.  The  expression  T.  u^H^CVjY)  +  Vc'  is  convex  in  V  >  0  and 


j=l  3  3 


Y  >  0. 


Condition  4  is  easily  shown  to  be  equivalent  to  a,  b,  and  c  providing 
n 

the  expression  7.  li.H^CVjY)  +  Vc'  is  quasi-convex.  The  fact  that 
3=1  3  3 

Condition  4  is  expendable  if  the  functions  H^(V,Y)  j=l,...,n  are 
convex  is  implied  by  d. 

Therefore,  our  new  set  of  conditions  imply  that  (V^Y0)  minimizes 

a  { 

~  u°K*(V,Y)  +  Vc'  subject  to  B(Y)  <  b,  V  >  0  and  Y  >  0,  i.e., 

i=l  J  J 


"  u<3H^(V°,YP)  +  V°c'  <  7  H^*CV,Y)  +  Vc' 

j=l  33  “  j=l  3  3 

for  all  V  >  0,  Y  >  0  and  B(Y)  <  b 


B(  Y°  )  1  b 

•which  is  Condition  3  of  the  previous  Section. 

Before  concluding  this  Section  some  observations  should  be  made 

n  Q  )  # 

about  the  properties  of  the  expression  Z  u^H^f(V,Y)  +  Vc'  in  the 

3=1  3  3 


27 


region  V  >  0,  Y  <:  0,  Condition  4  requires  that  this  expression  be 
quasi-convex,  i.e.,  every  set  R-  defined  as 

Z 


R-  =  {v,Y  |  2  U°H^(V,Y)  +  Vc'  <  z} 

*  j=l  3  J  ~  ’ 

must  be  convex.  If  the  functions  H^f(V,Y)  j=i,...,n  are  all  convex, 

then  by  virtue  of  the  fact  that  u°  ^  0  j=l,...,n,  the  expression 

n 

2  u.Htff(V,Y)  +  Vc'  must  also  be  convex.  However,  the  expression 
j=l  J  3 

is  not  necessarily  quasi-convex  if  the  functions  H|(V,Y)  ;=l,...,n 
are  each  quasi-convex.  To  establish  quasi-convexity  of 
n 

2  u-H,,f(V,Y)  +  Vc/  the  following  result  fr  m  Arrow  and  Enthoven  [1] 

j=l  3  J 

may  be  used.  A  function  F(X)  of  an  n-dimensional  vector  X  is  quasi7 

convex  if  D  <0  for  all  r=l,...,n  and  for  all  X  where  D  is  the 
r  =  r 

bordered  determinant : 


&F(X)/bX. 


.  . .  .£>F(X)/bX 


Dp  =  &F(X)/6X  b2F(X)/bX16X1  . . . .b2F(X)/5X1bX1 


j  bF(X)/bXr  ?>2F(X)/bXrbX1  ....b2F(X)/f,XriiXr  j  . 
If  it  is  possible  to  show  that  the  bordered  determinants  for  the 


function  F(X)  =  2  (i-H*f(V,Y)  +  Vc'  have  the  desired  sign  every- 

j=l  3  3 

where  in  the  region  V  >  0,  Y  >  0  then  our  Conditions  1  through  4  are 
sufficient  for  a  global  solution  to  the  constrained  min-max  problem. 
If  quasi-convexity  can  only  be  demonstrated  in  the  vicinity  of 
(V°,Y°),  then  a  local  solution  has  been  found  that  may  or  may  not 
also  be  a  global  solution. 
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V 


A  MODIFIED  DOUBLE  LAGRANGE  MULTIPLIER  METHOD 

In  general,  Lagrange  multiplier  methods  replace  constrained  problems 
with  unconstrained  problems  by  reversing  the  computational  roles  of 
resource  levels  and  imputed  resource  prices.  Such  methods  are 
’'fail-safe”  if  it  is  possible  to  show  that  a  solution  to  the  uncon¬ 
strained  problem  always  solves  a  constrained  problem  for  an  as¬ 
certainable  pair  of  resource  vectors.  The  value  of  the  method  as  a 
procedure  for  arranging  computations  is  entirely  attributable  to  the 
fact  that  explicit  constraints  can  frequently  be  an  expensive  compu¬ 
tational  nuisance. 

A  "fail-safe”  double  Lagrange  multiplier  method  for  constrained 
min-max  problems  with  mixed  strategies  that  conform  to  assumptions 
1  through  10  begins  with  an  arbitrary  choice  of  the  multiplier 
vectors  V  ^  0  and  W  >  0.  Next,  strategies  X°  j=l,...,n  for  the 

maximizing  player’s  tactics,  probabilities  j=l,...,n  and  a 
strategy  Y°  for  the  minimizing  player  are  found  such  that: 


n 

1.  X°  >  0  j=l . .  u°  >  0  j=l . .  E  =  1,  Y°  ^  0 

J  -  1  j=l  J 

2a.  For  every  tactic  j  and  each  element  X_.^  of  X^  : 

dH,(X°,Y°)  .. 

~bX^ - V  =  0  anS 

&H.(X°  Y°)  _ 

■  h^::-  - v  ° 0  if 


j* 


th 


where  a  ..  denotes  the  k  row  of  A  . . 
J*.  J 
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2b.  For  every  tactic  j: 


E  U°  [H,(X°  Y°)  -  V  A;x°']  >  H.(X°  Y°) 

3  3  -  3  3 


2  [H.(X°  Y°)  -  V  A'X°#]  =  H,(X°  Y°) 

j  =]_  ->  J  J  J  J  J  J 

For  each  element  Y,  of  Y: 

k 


n  rt  bH,(X°  Y°) 


V  0 

E  n. 
3=1  3 


T&b(y°)~|> 

L  5Yk  J  = 


0  and 


V  A'X°'  and 

J  3 

V  A'X°'  if  u°  >  0. 


=  0  if  y£  >  0 . 


The  resource  vectors  c  and  b  are  determined  after  the  fact  according 
to  the  formulae: 


c  =  £  U°X°A.  and  b  =  B(Y°) 
j=l  3  J  J 

and  the  solution  is  examined  to  verify  that  the  problem  is  not  de¬ 
generate  in  the  sense  that  Condition  4  is  not  satisfied.  It  should 

be  perfectly  apparent  that  the  strategies  X°  j=l,...,n,  probabilities 
o  o  3  _  - 

Uj  j=l,...,n,  the  strategy  Y  ,  and  the  vectors.  V  and  W  jointly  satisfy 

the  strong  set  of  sufficient  conditions  of  the  previous  Section. 

The  procedure  is  an  extrapolation  to  constrained  min-max  problems 
of  Everett’s  generalized  Lagrange  multiplier  method.  In  the  special 
case  of  a  single  tactic  the  procedure  reduces  to  the,  double  Lagrange 
multiplier  method  proposed  by  Pugh.  The  probabilities  then  drop 
from  view  and  Condition  2b  may  be  omitted. 

At  this  point  a  rough  intuitive  description  of  the  method  may 
be  helpful.  The  Lagrange  multiplier  vectors  V  and  W  bear  interpre¬ 
tation  as  prices  to  the  maximizing  and  minimizing  players  on  their 
resources .  These  vectors  provide  a  priori  rates  at  which  the  players 
may  convert  additional  units  of  resources  to  changes  in  the  expected 
value  of  the  payoff  function.  The  solution  value  of  the  Lagrange 
function : 
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2  U^H.(X°,V°)  -  V°  V  uVx°'  +W°  [B(Y°)]' 
j=l  J  3  3  -J=i  3  J  3 

also  has  an  economic  interpretation.  The  quantities  V°  2  q^A^X^' 

j=l  -1  3  J 

and  W°[B(V°)]'  are  the  opportunity  costs  in  payoff  units  of  the 
resources  used  by  the  maximizing  and  minimizing  players,  respectively. 
The  value  of  the  Lagrange  function  is  the  expected  payoff  to  the 
maximizing  player  modified  by  subtracting  the  opportunity  cost  of 
the  resources  he  must  expend  to  gain  this  payoff  and  by  adding  the 
value  of  the  defense  resources  used  to  limit  the  payoff. 

The  strategies  X°  j=i,...,n,  the  probabilities  u°  j=l,...,n, 
and  the  strategy  Y°  may  be  viewed  as  an  equilibrium  pair  in  a 
restricted  sense.  The  maximizing  player’s  mixed  strategy  is  optimal 
against  the  minimizing  player’s  strategy  but  the  reverse  relationship 
holds  in  only  a  partial  fashion.  The  strategy  Y°  is  not  necessarily 
strictly  optimal  against  the  specific  strategies  X°  mixed  with  probab¬ 
ilities  |i°.  Rather,  Y°  is  opcimal  against  a  specific  probabilistic 
mix  of  efficient  tactical  responses.  That  is,  the  fixed  probabilities 
d°  j=l,...,n  apply  to  strategies  X^  j=l,...,n  that  are  continuously 
altered  in  response  to  changes  in  the  strategy  choice  Y  and  the  imputed 
prices  of  the  minimizing  player’s  resources  V.  These  alterations 
have  the  effect  of  maintaining  the  strategies  X^  j=l,...,n  as 
economically  efficient  uses  of  resources  against  the  minimizing 
player’s  strategies.  The  specific  strategies  X^  j=l,...,n  are 
efficient  given  the  vector  V  and  the  solution  strategy  Y  . 
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VI 


ON  SUITABLE  PROBLEMS 

The  modified  double  Lagrange  multiplier  method  does  not  constitute 
a  uniformly  powerful  approach  to  all  constrained  min-max  problems 
with  the  requisite  characteristics.  Since  resource  levels  are 
treated  in  computations  as  an  output  rather  than  as  an  input,  the 
method  is  best  suited  to  problems  with  only  small  numbers  of  con¬ 
straints  and  to  investigations  of  solutions  over  a  range  of  resource 
levels.  Generally,  the  computations  required  to  yield  an  acceptable 
body  of  results  for  a  problem  will  tend  to  increase  geometrically  with 
the  number  of  constraints . 

The  method  is  principally  useful  as  a  means  for  addressing 
problems  that  may  be  partitioned  by  ’’cells”.  Let  h  =  l,...,m  be 
an  index  of  cells.  Each  cell  is  assumed  to  possess,  individually, 
all  of  the  components  of  a  constrained  min-max  problem  with  mixed 
strategies.  For  the  htn  cell,  then,  we  denote  these  components 
with  superscripts  as  follows: 

-  an  index  of  tactics  for  the  maximizing  player. 

-  a  strategy  for  the  maximizing  player’s  tactic  j. 

-  a  probability  of  use  of  Xj. 

-  a  strategy  for  the  minimizing  player. 

-  a  payoff  function  for  the  tactic  j . 

-  resource  use  by  the  maximizing  player  for  tactic  j. 

-  resource  use  by  the  minimizing  player. 

ie  maximizing  player  over  all  cells  is  a  specific 
combination  of  tactics  for  each  cell.  A  complete  list  of  such 
overall  tactics  would  consist  of  every  possible  combination  of  cell 
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tactics  and  can  be  quite  long  even  for  very  simple  problems.  Probabil¬ 
ities  for  these  overall  tactics  may  be  retrieved  from  the  probabilities 
assigned  to  the  cell  tactics  in  any  way  that  preserves  the  weights  iju 
as  conditional  probabilities.  For  example, 


h 

Hi  =  n  u!-  for 

3  h=l  3 


xj  =  {xj,xj . x”} 


corresponds  to  an  independent  random  drawing  of  the  cell  strategies 
1*1  h 

Xj  according  to  the  probabilities  .  An  overall  strategy  for  the 
minimizing  player  is  the  composite  strategy  Y  =  {y1, Y2, . . . ,Ym|  .  The 
expected  payoff  over  all  cells  is  the  sum  of  the  expected  payoff  in 
each  cell: 


The  modified  double  Lagrange  multiplier  method  makes  it  possiole 

to  solve  partitionable  problems  cell  by  cell.  That  is,  the  derivation 

,  o  ,  ,  o  , 

of  the  strategies  j=l,...,n  ,  the  probabilities  uu  j=l,...,n  , 

h° 

and  the  strategy  Y  for  any  one  cell  is  completely  separated  from 
the  derivation  for  any  other  cell.  This  fact  often  makes  it  possible 
to  solve  even  very  large  constrained  min-max  problems  by  partitioning 
them  into  cells  such  that  the  application  of  the  method  in  any  single 
cell  is  straightforward. 

Equally  important,  an  expectational  interpretation  of  the  payoff 
function  and  the  maximizing  player’s  resource  constraints  is  least 
likely  to  be  objectionable  with  problems  that  may  be  partitioned. 
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Suppose  that  the  maximizing  player's  solution  implies  an  independent 
random  drawing  of  strategies  in  each  cell.  This  random  process 
establishes  probability  distributions  for  the  actual  payoff  and  re¬ 
source  use  by  the  maximizing  player  that  would  result  from  a  play  of 
the  game  with  the  solution  strategies.  In  general  the  standard 
deviations  of  these  probability  distributions  will  tend  to  diminish 
relative  to  their  means  as  the  number  of  cells  is  increased.  This 
consideration  militates  strongly  in  favor  of  highly  separable  problems 
in  which  the  probable  variations  in  payoff  and  resource  use  become 
insignificant  in  relation  to  the  expected  values. 

The  modified  double  Lagrange  multiplier  method  is  a  general 
approach  to  formulating  and  solving  constrained  min-max  problems. 

It  is  not  a  computational  algorithm  in  the  same  sense  as  the  simplex 
algorithm  or  Newton’s  method.  In  deriving  solutions  that  satisfy  the 
Ponitions  of  the  previous  section,  the  analyst  is  left  to  his 
own  devices.  The  method  applies  equally  well  with  clever  mathe¬ 
matical  derivations  and  unsubtle  searches  on  a  digital  computer. 

The  choice  of  a  computational  route  for  any  particular  problem 
must  be  made  on  the  basis  of  the  specific  characteristics  of  the 
problem.  Moreover,  we  have  provided  no  grounds  to  warrant  the 
belief  that  solutions  to  all  suitable  constrained  min-max  problems 
can  be  found  by  applying  the  method.  The  conditions  we  have  derived 
are  sufficient  but  they  may  not  all  be  necessary. 
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VII 


AN  ABM/SHELTER  DEPLOYMENT  PROBLEM 

A  simple  ABM/shelter  deployment  problem  can  be  solved  as  an  illustra¬ 
tion  or  the  essential  features  of  the  modified  double  Lagrange  multi¬ 
plier  method.  The  game  pits  warheads  against  terminal  interceptors 
and  blast  shelters  allocated  to  cities.  The  maximizing  player’s 
strategies  are  targeting  plans  in  which  each  of  the  cities  is  assigned 
a  specific  number  of  warheads.  Strategies  for  the  minimizing  player 
consist  of  assignments  of  interceptors  and  shelters  to  cities.  The 
interceptors  are  taken  to  be  perfect  in  the  sense  that  any  warhead 
engaged  by  the  defense  is  certain  to  be  destroyed.  We  shall  also 
assume  that  each  shelter  provides  one  resident  of  the  city  with  com¬ 
plete  protection  from  all  weapons  effects. 

The  payoff  of  the  game  is  the  total  number  of  fatalities  anti¬ 
cipated  from  an  attack.  In  any  single  city  fatalities  occur  among 
the  unsheltered  population  and  are  computed  as  an  exponential 
function  of  the  number  of  detonating  v/arheads.  Both  players  are 
constrained  by  a  single  scarce  resource.  In  the  case  of  the  maximizing 
player,  his  arsenal  is  presumed  to  consist  solely  of  ballistic 
missiles  capable  of  delivering  only  a  limited  number  of  warheads. 

The  defense  is  restricted  to  deployments  that  do  not  require  expendi¬ 
tures  in  excess  of  a  given  budget.  His  cost  function  is  linear. 

Since  the  maximizing  player  observes  the  defense  before  choosing  a 
targeting  plan  we  have  a  constrained  min-max  problem.  (This  formu¬ 
lation  of  the  ABM/shelter  deployment  problem  is  not  intended  to  be 
more  than  illustrative.  An  application  of  the  method  to  a  more 
sophisticated  and  realistic  version  of  the  problem  has  been  made 
elsewhere  by  the  author. ) 

The  problem  may  be  partitioned  by  cities.  Let  h=l,...,m  be  an 
index  of  cities.  For  the  htn  c.ty  the  maximizing  player  exercises 
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one  of  two  tactical  options.  He  may  choose  not  to  attack  the  city 
at  all.  This  we  shall  call  the  null  attack.  The  second  tactic  may 
be  labeled  an  exhaustion  attack.  In  an  efficient  use  of  this  tactic 
the  defense’s  interceptors  are  exhausted  and  the  maximizing  player 
continues  to  target  warheads  until  the  last  warhead  is  just  worth  the 
damage  it  will  cause.  Mathematically,  the  components  of  the  game  for 
the  htn  city  are: 

Xh  -  number  of  warheads  for  the  exhaustion  attack, 

li*1  -  probability  of  using  the  exhaustion  attack, 

Yn  -  a  two-element  vector  as  follows : 

-  the  number  of  terminal  interceptors 


Y^  -  the  number  of  blast  shelters, 

Ph  -  the  population  of  the  city, 

sh  -  a  damage  function  parameter, 
a  -  cost  per  interceptor,  and 
0  -  cost  per  shelter. 


The  payoff  function  for  the  null  attack  is : 

hJco^)  =  0 

and  for  the  exhaustion  attack: 


HjCxb.yb) 


J° 

J(Ph  -  yS)Ci  - 


exp  {-  sh(Xh  -  yJ)}] 


if  Ph  -  Y^  <  0 
if  ph  -  y£  >  c . 
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These  two  functions  are  graphed  in  the  figure  below  for  P^1  -  >  0 

and  >  0: 


We  have  no  vector  for  the  null  attack,  so: 

Hjcv.yb)  =  o. 

For  the  exhaustion  attack  we  have: 


H„(VYI..J°  ifPh-V*<0 

2  *  (  Supxh  |(Ph  -  Y^)[l  -{exp  -sh(Xh  -  yJ  } ]  -  VXh|xh  *  oj 

if  Ph  -  Y^  >  0 

If  V  =  0,  the  supremum  takes  the  value  (Ph  -  y!?);  if  V  >  0,  the 
supremum  is  a  maximum.  Differentiating  with  respect  to  X  yields: 

6H5(Xh,Yh)  -  VXh 


=  sh(Ph  -  y£)  exp  {-sh(Xh  -  yJ)|  -  V. 


At  a  maximum: 

sh(Ph  -  Y$)  exp  {  -sh(Xh  -  yJ)}  -  V  1  0  and 
sh(Ph  -  yJ)  exp  {  -sh(Xh  -  Yj)}  -  V  *  0  if  Xh  >  0. 
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Therefore  ? 


0  if  Ph  -  %  0 

HJCV,**1)  =  (Ph  -  v5)[l  -  exp^}  ]  if  0  <  Ph-Y$  ^exp^sN^/s*1 
(Ph  -  Yjj-VCl+logjs^-^J/V^/s^VYj 

if  ph  “  *2  ^  Vexpl-s^j/sb. 

The  expected  payoff  in  the  nation  as  a  whole  is  the  expression: 

0  if  Ph  -  y£  <  0 

J  h  2  = 

E  M 

(Ph  -  y£)[1  -  exp|-sh(Xh  -  yJ)J]  if  Ph  -  Y*  >  0. 


The  maximizing  player’s  expected  use  of  warheads  is: 

£  uV 

h=l 

and  the  cost  of  the  ABM/shelter  deployment  is  the  linear  expression: 

i  [ay!!1  +  by!)]. 

h=l  1  i 

It  may  be  readily  verified  that  the  problem  has  all  of  the  character¬ 
istics  (Assumptions  1  through  10)  that  are  necessary  for  an  appropriate 
application  of  the  modified  double  Lagrange  multiplier  method. 

To  apply  the  modified  double  Lagrange  multiplier  method  we  select 
multipliers  9  >  0  and  #  >  0  and  then  choose  strategics  and  probabili¬ 
ties  that  constitute  an  equilibrium  pair  in  the  restricted  sense  of 
the  previous  sections.  Omitting  the  superscript  h,  Cor.ni^iwns  1,  2a, 
2b,  and  3a  reduce  to  the  following  system  of  conditional  equations  and 
inequalities  for  each  city: 


X°  =  0  if  (P-Y^)-9/s  1  0 

X°  =  log|s(P-V^)/vJ/s  +  V®  if  (P-Y^)-V/s  >  0 

u°  =  0  if  (P-Y°)  -  tf[l  +  log  js(P-Y^)/$  j]/'s  -  9Y®  <  0 
0  ^  U°  1  1  if  (P-Y§)-fl[l+log{s(P-Y§)/0}]/s  -  9y£  =  0 
u°  -  1  if  (P-Y^)-9[l+log js(P-Y^)/9j]/s  -  9yJ  >  0 
Y^  =  0  if  u°  [-V]  +  ft  a  >  0 

Y°  \  0  if  u°  C-V]  +  W  a  =  0 

Y°  =  0  if  u°  [V/s(P-Y^)  -  1]  +  W0  >  0 

Y°  ^  0  if  d°  [9/s(P-Y^)  -  1]  +  W8  =  0. 

There  are  five  distinct  types  of  solutions  to  this  system  dis¬ 
cussed  below  as  cases  1  through  5. 

Case  1:  sP  <  V 


|i°Hr(X°,Y°)  =  0 

u°X°  =  0 

ciY^  +  8Y^  =  C. 

The  undefended  city  is  not  an  attractive  target  to  the  maximizing 
player.  Therefore,  there  is  no  reason  to  install  interceptors  or 
build  shelters.  In  all  of  the  remaining  cases,  the  undefended  city 
presents  a  worthwhile  target  to  the  maximizing  player. 
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Case  2:  sP  ^  9,  Wa  >  V,  wp  >  1  -  9/sP 

X°  =  log|sP/v|/s 

Ii°  =  1 

vj  =  0 

=  0 

u°h2(x°,y°)  =  P  -  9/s 
d°X°  =  log  jsP/V  j/s 
aY°  +  BY^  =  0. 

Although  the  city  is  subject  to  attack,  the  minimizing  player  finds 
that  it  is  best  to  leave  the  city  undefended.  The  reasons  for  this 
are,  first,  that  the  opportunity  cost  in  lives  of  an  interceptor, 

Qa,  exceeds  the  value  in  lives  of  the  single  warhead  it  destroys,  V. 

Second,  the  opportunity  cost  in  lives  of  a  shelter,  Wp,  exceeds  the 
expected  saving  in  lives  of  installing  the  shelter,  1  -  V/sP. 

Case  3:  sP  ^  V,  Wa  >  V,  Wp  ^  1  -  V/sP 

X°  =  log jl/(l  -  ftp)} 

o  _  . 
d  =  1 

y£  =  0 

y£  =  P  -  v/s(l  -  wp) 
d°H2(X°,Y°)  =  V  Qp/s(l  -  Qp) 
d°X°  =  log|l/(l  -  ftp)} 
aY^  +  BY^  =  BP  -  9p/s(l  -08). 
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In  this  case  the  blast  shelters,  but  not  the  interceptors,  are  worth 
the  cost.  The  optimum  number  of  shelters  minimizes  the  Lagrangian: 

U°Hf(V,Y)  +  W0Y2. 

Differentiating  at  Y  =  Y°  we  find  that: 

U°[  -  1  +  9/s(P  -  Y^)3  +  We  =  0. 

At  the  margin  no  increase  or  decrease  in  the  number  of  shelters  leads 
ultimately  to  a  decrease  in  the  overall  payoff.  Since  the  offense  may 
move  warheads  among  cities  to  maintain  a  marginal  return  of  9  per  war¬ 
head,  the  net  number  of  lives  saved  by  an  additional  shelter  in  the 
city  is  [1  -  9/s(P  -  Y^)].  This  just  balances  the  opportunity  cost  in 
lives  of  the  shelter,  fl0. 

Case  4:  sP  ^  9,  fla  %  9,  553  >  (Qa/9)(1  -  9/sP) 

X°  =  P/9  -  1/s 
H°  =  Wa/9 

Y°  =  P/9  -  [1  +  log jsP/9  |]/s 


U0H2(X0,Y°)  =  (fla/9)[P  -  9/s] 

d°X°  =  (Wa/9) [p/9  -  i/s] 

aY°  +  0Y2  ~  aP/V  -  a[l  +  logjsP/v|]/s, 

The  city  is  allocated  interceptors  but  not  shelters  since  the 
opportunity  cost  of  a  shelter,  flg,  exceeds  the  anticipated  saving 
in  lives  from  an  attack,  1  -  9/sP,  multiplied  by  the  probability 
of  an  attack,  fla/9.  The  optimum  number  of  interceptors,  Y^,  is 
determined  not  by  a  margin  calculation  as  in  case  3,  but  is  set 
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instead  at  a  level  that  makes  attacking  the  city  at  all  an  indif¬ 
ferent  proposition  to  the  maximizing  player.  That  is: 

H*(V,Y°)  =  HgCV, Y°)  =  0. 

The  null  attack  and  the  exhaustion  attack  tactics  are  probabilistically 
mixed  with  probabilities  (1  -  VJa/V)  and  $a/V  such  that  neither  an 
increase  nor  a  decrease  in  the  number  of  interceptors  will  reduce 
the  value  of  the  Lagrangian: 

(1  -  u°)H*(V,Y°)  +  U°H*(V,Y°)  +  WaY^  . 

Differentiating  with  respect  to  Y^  and  evaluating  the  derivative  at 
Y°  we  have : 


o 

d 


&h*(9,y°) 


— [-  V]  +  Wa  =  0. 

.9  / 


A  rough  intuitive  explanation  for  the  choice  of  u°  =  V?a/V  may  be 
provided  by  considering  the  role  of  the  probability  p°  in  ascertaining 
the  resource  level  c.  If  we  choose  d  <  Wa/V,  then  the  resulting 
arsenal  is  nor  sufficiently  large  to  discourage  a  defender  from 
decreasing  the  number  of  interceptors  in  every  defended  city.  On 
the  other  hand  d  >  Wa/V  yields  a  total  number  of  warheads  so  large 
that  an  increase  in  the  number  of  interceptors  at  every  defended  city 
must  produce  a  return  in  the  imputed  value  of  intercepted  warheads 
that  on  the  average  exceeds  Wa  per  additional  interceptor. 

Case  5:  sP  >  V,  Wa  <  V,  W0  -g  (Wa/V)(l  -  V/sP) 


X°  =  a/s(a  -  V8)  -  1/s 


U°  =  Wa/V 

Yj  =  a/s(a  -  9g)  -  1/s  -  log{a/(a  -  9s)|/s 
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Y°  =  P  "  Va/s(a  -  V8) 

U°H2(X°,YD)  =  (Wa/V)[Va/s(a  -  Vg)  -  V/s] 

U°X°  =  (fla/9)[a/s(a  -  Vg)  -  1/s] 

aY°  +  gY^  =  8P  -  a  log{a/(a  -  Vp)J/s. 

This  last  arid  most  interesting  type  of  solution  is  a  mixture  of  active 
and  passive  defense  systems.  An  heuristic  explanation  of  the 
solution  is  a  straightforward  composite  of  the  descriptions  provided 
for  cases  3  and  4. 

The  various  values  of  the  multipliers  and  their  association  with 
the  five  types  of  solutions  are  depicted  in  the  figure  below. 


tlQ  t"lO 

Once  a  number  of  warheads,  X  ,  probability,  u  ,  number  of 

interceptors,  Y^°  ,  and  number  of  shelters,  Y^°,  have  been  determined 

for  every  city,  the  last  step  of  the  modified  double  Lagrange  multi¬ 
plier  method  is  a  simple  summation  of  expected  payoffs,  expected  use 
of  warheads,  and  defense  expenditures.  In  the  normal  course  of  an 
operations  research  exercise  this  would  be  repeated  for  many  combina¬ 
tions  of  values  for  the  multipliers  V  and  W. 
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Unfortunately,  it  is  only  possible  to  show  that  the  strategies 
derived  using  the  modified  double  Lagrange  multiplier  method  are  local 
solutions  to  the  ABM/shelter  problem.  An  examination  of  the  bordered 

n 

determinants  for  the  expression  F(X)  =  I  ia.H1i?(V,Y)  +  Vc'  indicates 

j=l  :  3 

they  will  all  be  negative  wherever: 

||  =  Z  logjsh(Ph  -  Yfj})/vJ/sh  -  vj  +  Xh0  0. 

This  derivative  is  equal  to  zero  at  the  solution  point  V  =  V,  Y  =  Y° 
but  is  not  generally  nonpositive  over  the  entire  region  V  ^  0,  Y  ^  0. 
It  can  be  easily  shown  that  part  c  of  Condition  4  is  met  if  at  least 
one  city  is  both  targeted  and  defended. 
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VIII 


CONCLUDING  REMARKS 

There  are  three  major  criticisms  that  can  be  made  of  the  approach 
to  constrained  min-max  problems  offered  in  this  paper: 

1.  The  wrong  problem  has  been  solved.  A  solution  to  a  con¬ 
strained  min-max  problem  that  implies  a  mixed  strategy  for  the 
maximizing  player  can  have  no  direct  application. 

2.  Ten  strong  assumptions  are  needed  to  obtain  tne  set  of 
sufficient  conditions  that  finally  provide  the  basis  for  the  modified 
double  Lagrange  multiplier  method.  Also,  quasi-convexity  of  the 

n 

expression  Z  (j.°H':(V,Y)  +  Vc'  (Condition  4)  is  not  usually  easy 
j=l  J  J 

to  establish, especially  over  the  entire  region  V  :>  0,  Y  ^  0. 

3.  As  a  practical  matter,  the  use  of  the  method  is  limited  to 
partitionable  problems  with  few  constraints  and  to  crude  explorations 
of  solutions  over  a  range  of  resource  levels . 

My  response  to  these  criticisms  is: 

1.  The  right  problem  has  proved  to  be  quite  difficult  to 
address  on  a  practical  level  with  existing  computational  algorithms. 
The  introduction  of  mixed  strategies  and  the  redefinition  of  an 
equilibrium  pair  appears  to  me  to  be  the  least  damaging  alteration 
that  leaves  a  readily  solveable  version  of  the  constrained  min-max 
problem.  Solutions  to  this  wrong  problem  will  often  provide  close 
approximations  or  at  least  bounds  for  solutions  to  the  right  problem. 

2.  The  ten  assumptions  we  have  used  are  probably  not  the  only 
assumptions  on  which  a  Lagrange  multiplier  method  for  constrained 
min-max  problems  can  be  based.  For  example,  a  fail-safe  version  of 
the  method  could  proceed  directly  from  the  first  set  of  sufficient 
conditions  providing  Condition  3a  of  Section  III  can  be  weakened  to: 
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n 

Z  u°H^(V°  ,  Y°)  +  W°[B(Y°)]* 
j=l  J  3 

n 

1  £  d°H.(V°,Y)  +  W°[B(Y)]'  for  all  Y6T. 
j=l  3  3 

I  doubt  that  there  are  many  problems  for  which  this  cannot  be  done. 
Also,  it  may  be  possible  in  many  circumstances  to  define  the 
functions  j=l,...,n  more  generally  as: 

Hj(V, Y)  =  Supx  I^CX^Y)  -  V[Cj(Xj)]'  |  Xj6Sj|  . 

Further  research  may  also  reveal  that  Condition  4  of  Section  IV  may 
be  weakened  particularly  with  respect  to  the  quasi-convexity  of 
n 

£  U^HvCVjY)  +  Vc*  over  the  entire  region  V  >  0,  Y  >  0. 
j=l  3  3  " 

3.  The  practical  limits  on  the  usefulness  of  our  method  are 
about  the  same  as  those  that  apply  to  Everett’s  original  generalized 
Lagrange  multiplier  method.  Many  problems  of  interest  are  net 
well-suited  for  applications  of  Lagrange  multiplier  methods.  How¬ 
ever,  some  of  the  problems  that  are  well-suited  are  just  about 
impossible  to  address  in  any  other  way.  For  example,  the  author 
has  been  engaged  in  a  study  of  mixes  of  area  ABM,  terminal  ABM, 
and  blast  shelters  to  defend  the  nation's  urban  population  from 
strategic  missile  attacks.  The  admissible  offense  strategies  involve 
a  variety  of  weapons  and  decoys  deployed  to  leak  through,  suppress, 
exhaust,  and  detonate  above  the  terminal  defenses.  Despite  its 
apparent  complexity  the  problem  can  Le  formulated  as  a  min-max 
problem  with  constraints  on  the  offense  use  of  throw  ’weight  and  the 
defense  use  of  budget  funds.  The  problem  can  be  partitioned  by 
cities  and  rapidly  and  repeatedly  solved  for  different  offense  and 
defense  resource  levels  simply  by  varying  two  multipliers.  For 
entirely  practical  reasons  no  other  approach  to  this  constrained 
min-max  prob.lem  is  currently  feasible. 
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